The four jet production at LHC and Tevatron in QCD. 
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We demonstrate that in the back-to-back kinematics the production of four jets in the collision 
of two partons is suppressed in the leading log approximation of pQCD, compared to the hard 
processes involving the collision of four partons. We derive the basic equation for four-jet production 
in QCD in terms of the convolution of generalized two-parton distributions of colliding hadrons in 
the momentum space representation. Our derivation leads to geometrical approach in the impact 
parameter space close to that suggested within the parton model and used before to describe the 
four-jet production. We develop the independent parton approximation to the light-cone wave 
function of the proton. Comparison with the CDF and DO data shows that the independent parton 
approximation to the light-cone wave function of the proton is insufficient to explain the data. We 
argue that the data indicate the presence of significant multiparton correlations in the light-cone 
wave functions of colliding protons. 
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In spite of extensive theoretical and experimental work, 
various aspects of the high-energy hadronic collisions at 
the Tevatron and LHC arc still poorly understood. This 
is especially true for the multijet production which is of 
a paramount importance for the understanding of pQCD 
dynamics at high-energy colliders, and for the search of 
new particles. The topic of multiparton interactions is 
now one of the focal points of studying pQCD and build- 
ing an adequate basis for modeling the final states at the 
LHC. In particular, description of multiparton interac- 
tions requires treatment of a significant inbalance of the 
momenta of the jets (presence of the Sudakov form fac- 
tors). In this paper we summarize the first steps of the 
program to address these issues. Among the original re- 
sults of the paper are the derivation of the formulas in 
the leading logarithmic approximation for production of 
4 jets. Our key finding is that it is possible to isolate the 
kinematics where the leading twist processes 2->4 are 
not enhanced. This result will allow one to improve the 
reliability of the Tevatron studies of the four-jet produc- 
tion in the multiparton kinematics and point out direc- 
tions for the corresponding analysis at the LHC. 

Another critical issue is the formulation of the prob- 
lem in terms of the generalized two-parton distributions 
in the momentum space representation and introduction 
of the mean field approximation for this object. This 
new formulation is well suited for the more detailed stud- 
ies which are now under way. In addition it establishes a 
link with the original formulation in the coordinate space 
[Tj-TToj , and resolves an issue of the value of the strength 
of the double interaction within this approximation. Pre- 
viously there was a question whether a conclusion of Ref . 
[H , that the observed rate is a factor of 2 larger than the 
theoretical prediction, can be due to uncertainties related 



to the many Fourier transforms required to convert the 
HERA data to the experimental number. A new formula- 
tion, though mathematically equivalent, has completely 
resolved this issue. This poses serious constraints on the 
Monte Carlo models of pp scattering at collider energies 
which are not satisfied by many of the current models. 

These issues are of broad interest, both theoretical and 
experimental. 

The standard approach to the multijet production is 
the QCD improved parton model. It is based on the as- 
sumption that the cross section of a hard hadron-hadron 
interaction is calculable in terms of the convolution of 
parton distributions within colliding hadrons with the 
cross section of a hard two-parton collision. An applica- 
tion of this approach to the processes with production of 
four jets implies that all jets in the event are produced 
in a hard collision of two initial state partons. 

The recent data of the CDF and DO Collaborations 
pU4l3j do not contradict to the dominance of this mech- 
anism in the well-defined part of the phase space. At 
the same time these data provide the evidence that there 
exists a kinematical domain where a more complicated 
mechanism becomes important, namely the double hard 
interaction of two partons in one hadron with two partons 
in the second hadron. 

Within the parton model picture, the four jets pro- 
duced this way should pair into two groups such that the 
transverse momenta of two jets in each pair compensate 
each other. In what follows we refer to this kinematics 
as back-to-back dijet production. We consider the dijets 
for the case 



sf 3 = tiu + htY « fu - jit, « At - Jit, (i) 

where S is the total transverse momentum of the dijet and 
jit the transverse momentum of an individual jet (see 
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Fig. [I}. The hardness condition S 2 3> R~ 2 is implied, 
with R the characteristic hadron size (nonperturbativc 
scale). The events with inbalances S 2 < R~ 2 give a small 
contribution both to total and differential cross sections, 
since they are suppressed by the Sudakov form factors. 
For a detailed discussion of this issue see review [13] • 

Importantly, in this kinematical region the hard scat- 
tering of four partons from the wave functions of the col- 
liding hadrons remains the dominant source for four-jet 
production even when the pQCD parton multiplication 
phenomena are taken into account. 

The reason for that is the following. When the two 
partons from each hadron emerge from the initial state 
parton cascades and then engage into double hard scat- 
tering, the resulting differential distribution of the fi- 
nal state jets lacks the double back-to-back enhancement 
factor da oc S^ 3 6^4 which is there in the case of two 
independent hard scatterings [Til ]. For the two-parton 
scattering, the characteristic perturbative enhancement 
da oc 6~ 2 results from a coherent enhancement of the 
amplitude due to integration over a large transverse disk, 
p 2 ~ 5~ 2 3> ii" 1 . The two partons that originate from a 
perturbative splitting form a relatively compact system 
in the impact parameter space, so that the double hard 
interaction of such pairs produces only a single pertur- 
bative enhancement factor, {6\ 3 + 524)~ 2 , which does not 
favor the back-to- back dijet kinematics (fTJ). The distri- 
bution of four jets so produced is much more isotropic 
and can be suppressed by choosing proper kinematical 
cuts. 

So, the aim of this letter is to consider the four-jet pro- 
duction in the hard collisions of four initial state partons. 
We show that the cross section of back-to-back dijet pro- 
duction is calculable in terms of new nonperturbative ob- 
jects — the generalized two parton distributions (2GPD) 
The properties of the 2GPD can be rigorously studied 
within QCD. In particular, we report here the derivation 
of the geometric picture for multiple parton collisions in 
the impact parameter space. 

In the kinematical domain (TTJ the direct calculation of 
the light-cone Feynman diagrams (momenta of the par- 
tons in the initial and final states are shown in Fig. [1]) 
using the separation of hard and soft scales shows that 
the four — > four cross section for the collisions of hadrons 
"a" and "b" has the form: 




FIG. 1: Kinematics of double hard collision - momenta of the 
colliding partons in \in > and < out\ states. 



contribution.). Remember that the light-cone fractions 
Xi are actually fixed by the final jet parameters and the 
energy momentum constraints. 

With account of the radiative pQCD effects, in full 
analogy with the "DDT formula" for two-body collisions, 
the differential distribution ^ acquires Sudakov form 
factors [3 HI] depending on the logarithms of the large 
ratios of scales, j 2 /S 2 , and the 2 GPDs become scale de- 
(5| 4 . It should be mentioned that 



pendent: pi ~ 5 2 3l p\ 
the structure of the final formula depends on what one 
actually measures in the experiment — energetic single 
particles with large transverse momenta in the final state 
or "jets" — and on how the jets are precisely defined. A 
more detailed account of the pQCD effects will be given 
in a future publication [15j . 

For brevity we will not write explicitly the virtuality 
scales of the 2GPD and will use the form: D(xi,X2, A). 
Note that these distributions depend on the new trans- 
verse vector A that is equal to the difference of the mo- 
menta of partons from the wave function of the colliding 
hadron in the amplitude and the amplitude conjugated. 
Such dependence arises because the difference of parton 
transverse momenta within the parton pair is not con- 
served. The integration limits in Xi, t are subject to stan- 
dard limits determined by experimental kinematic cuts. 

Within the parton model approximation the cross sec- 
tion has the form: 



a i (x 1 ,X2,X3,x 4: ) = J 



d 2 A 

(2tt) 2 



D a (x 1 ,x 2 ,pl,pl, A) 



xD b (x 3 ,x 4 ,p 2 1 ,p 2 ., 



-1 



da 13 da 2i - 

x — - —dtidt 2 . 

dt\ dti 



(2) 



Here D a (xi, X2,p 2 ,P2, A) are the new 2 GPDs for 
hadrons "a" and "b" defined below. (In the following we 
will consider the case of pp collisions and omit the sub- 
scripts a and b. Summing over collisions of various types 
of partons is implied. In practice however we will keep 
hard scattering of gluons only since it gives the dominant 



(74 = a\a2l^R\ 



(3) 



where a\ and 02 are the cross sections of two independent 
hard binary parton interactions. The factor 7ri?? nt which 
in principle depends on Xi characterizes the transverse 
area occupied by the partons participating in the hard 
collision. (In the experimental [ill, 113 and some of the 
theoretical papers this factor was denoted as an effective 
cross section. Our Eq.|3]below shows that such wording is 
not satisfactory since 7ri?? nt does not have the meaning of 
the interaction cross section.) The data [Ill4l3j indicate 
that 7ri2p nt is practically constant in the kinematical range 
studied at the Tevatron. 
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Eq. [2] leads to the general model independent expres- 
sion for 



1 



d 2 l 



D(xx,x 2 , 



(2tt) 2 D(x 1 )D(x 2 )D(x 3 )D{x i ) 



(4) 



in terms of the 2GPDS. Here D{xj) are the corresponding 
structure functions. Note that hereafter we do not write 
dependence of 04 and R 2 nt on the light-cone fractions Xi 
explicitly. 

The 2GPDS are expressed through the light-cone wave 
functions of the colliding hadrons as follows. Suppose 
that in a four — > four process the two partons in the nu- 
cleon in the initial state wave function have the transverse 
momenta k±,k 2 . Then in the conjugated wave function 
they will have the momenta k\ + A,k 2 — A . This is 
because only the sum of parton transverse momenta but 
not the difference is conserved. 

The relevant 2GPDS are: 



0(pI - k\ 



n 

i#l,2 



d 2 h 

W) 2 



x ip n (xi,ki,X2,k2,.,ki,Xi..) 



— y — y — y — y -* 
x ip+(xt,ki + A,x 2 ,k 2 - A,x 3 ,k 3 ,...) 



(5) 



i=l 



Note that this distribution is diagonal in the space of all 
partons except the two partons involved in the collision. 
Here tp is the parton wave function normalized to one 
in a usual way. An appropriate summation over color 
and Lorentz indices is implied. In the case of kinematics 
1 3> X\ > x 2 we expect only distributions without the 
spin flip to be important. 

Let us stress that it follows from the above formulas 
that in the impact parameter space these GPDs have a 
probabilistic interpretation. In particular they are posi- 
tively definite in the impact parameter space, cf. Eq. 1111 
Note that in the same way one can introduce the N- 
particle GPD, Gjv, which can be probed in the produc- 
tion of N pairs of jets. In this case the first N arguments 
ki in Eq. [5] are shifted by A* subject to the constraint 
J2i A, = 0. So the cross section is proportional to 



cr 2 jv oc 



5 (2*) 



' £> (Ai,...Ajv) 



D b (t u 



i=N 



(6) 



These GPDs can be easily rewritten in the form of the 
matrix elements of the operator product. For example: 



13(A) = <N\J d^x^x^x-i 

x G? + (x 1 )G b J+ (x 2 )G? + (x 3 )G b J+ (x A ) 

x exp(ipj(x\ — x 3 )~ + ip 2 {x 2 - £4)" 

+ iA t (x 4 -x 3 )t)\N >, 



(7) 



calculated at the virtualities p\ , p\ at fixed A . Here we 
gave an example for the most relevant case of gluons with- 
out a flip in color and spin spaces. In general a number 
of distributions can be written, depending on different 
contractions of transverse Lorentz indices and color in- 
dices. The classification of the relevant distributions is 
the same as the classification of the quasipartonic opera- 
tors in Ref. [17| . Note that the presence of the transverse 
external parameter A does not change the classification, 
since the corresponding new structures will be strongly 
suppressed at high energies. We wrote the operator ex- 
pression in the light-cone gauge. In an arbitrary gauge 
we shall need Wilson lines W(C) connecting points with 
contracted color indices. 

In the approximation of uncorrelated partons it follows 
from Eq. \5\ that 

D{ Xi ,x 2 ,pIpI A) = G( Xl ,pl A)G(x 2 ,pl A), (8) 

where G(x, A) are conventional one-particle GPDs. 
These GPDs can be approximated as Gn{x, Q 2 , A) = 
Gn{x 1 Q 2 )F 29 {A), where F 2g (A) is the two-gluon form 
factor of the nucleon extracted from hard exclusive vec- 
tor meson production (we suppress here the dependence 
of F 2g on x) [l8| and Gn{x, Q 2 ) conventional parton dis- 
tribution of a nucleon. (Here Q 2 is the virtuality due to 
the radiation, cf. discussion after Eq. [2j) Thus : 



irRi 



d 2 A 



a 

28tt' 



(9) 



Here at the last step we used the dipole fit F 2g (A) = 
l/(A 2 /m 2 + l) 2 to the two-gluon form factor (m 2 (x ~ 

0.03, Q 2 = 3GeV 2 ) « l.lGcV 2 ). Using the transverse 
gluon radius of the nucleon we obtain 



7/2r 



dF 2g (t)/dt t=0 . (10) 



This result coincides with the one for the area 7ri? 2 nt ob- 
tained earlier in [J] using the geometric picture in the 
impact parameter space. That derivation involved taking 
the Fourier transform of the two-gluon form factor and 
calculating a rather complicated six-dimensional integral 
which could potentially lead to large numerical uncertain- 
ties. The form of Eq. [TU] clearly indicates that numerical 
uncertainties are small. 
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It was emphasized in [4J that the experiments on four- 
jet production report a smaller value of 7ri? 2 nt as com- 
pared to the one obtained above in the independent par- 
ticle approximation (though the issue of how well the 
contribution of the 2 — » 4 processes was subtracted still 
remains, cf. discussion in beginning of the paper). It is 
at least a factor of 2 smaller — that is a four-jet cross 
section is a factor of 2 larger — than Eg. [TU1 gives. (The 
GPDs for sea quarks appear to decrease with A some- 
what faster, resulting in a smaller l/7i\R 2 nt , see discussion 
in[l9|.) 

It follows from Eq. 0] that the value of Rf nt is deter- 
mined by the range of integration over A. Hence the 
characteristic A in the integral measures the effective dis- 
tance between the parton pairs (which in principle may 
differ for different flavor combinations). According to 
the above evaluation within the independent parton ap- 
proximation the integral for 1 / R? nt is dominated by small 
A 2 ~ 0.1 to 2 ,. The contribution of large A is suppressed 
by the two-gluon form factor of a nucleon. This reasoning 
indicates the important role of interparton correlations. 
In other words, the integral over A is effectively cut off 
by a scale of the nonperturbative correlations. Such cor- 
relations naturally arise in nonperturbative QCD regime 
in a number of nucleon models, such as constituent quark 
model (gluon cloud around constituent quark) or 
string model (gluon structure of string) |2(|. The de- 
tailed analysis of the additional correlations due to the 
hard- soft interplay will be reported elsewhere (la ]. 

Let us now show that results obtained in the paper 
lead to the geometric picture in the impact parameter 
space mentioned above [H-[Ic||. 

The first step is to make transformation into coor- 
dinate space i.e., to make the Fourier transform from 
variables ki in Eq. [5] to coordinates pi. Performing in- 
tegration over ki we obtain that transverse coordinates 
of partons in the amplitude and the amplitude conju- 
gated are equal pi = pf. In the calculation we use 
the fact that upper limit of integration over k 2 is very 
large compared with the inverse hadron size. The next 
step is to perform integration over A which produces 
Hpi - fc- fc + Pi) = J d 2 B6(p 1 - p 3 ~ B)S{p 2 - Pa- B). 

The delta functions express the fact that within the ac- 
curacy 1/pt where pt is the hard scale, the interactions of 
partons from different nucleons occur at the same point. 
B is the relative impact parameter of two nucleons. 



The expression for the cross section in the impact pa- 
rameter space has the form which corresponds to the ge- 



ometry of Figj2] 
r 

<j 4 = d 2 Bd 2 p 1 d 2 p 2 d 2 p 3 d 2 p4 : D(x 1 ,x 2 ,p 1 . / p 2 ) 

x D{x 3n x i ,p 3 ,p i )8{pi - (B + p 3 ))S(p 2 - (B + pi)) = 
= J d 2 Bd 2 p 1 d 2 p 2 D(x 1 ,x 2 ,pi,P2) 

x D(x 3 ,x 4 ,-B + p u -B + p 2 ). (11) 

Here the 2GPD in the impact parameter space represen- 
tation is given by 

D (x\,X2,p\,p 2 ) = 
n—oc „ i—n 

= ^2 j Y\_[dxid 2 pi]i) n {xi,pi 1 x 2 ,p 2 ,...x ll p l ,) 

n=3 J i>3 

i—n 

X ^i(xi,Pi,X2,P2,-,Xi ) pi ) ...)6(^2x i pi). (12) 

i=l 

where the delta function expresses the center of mass 
constraint X);=i x iPi = 0- This is analogous to the 
case of single parton GPDs, see [2ljj . The functions 
i]j{xi, pi,x%, p2, ■■■) are just the Fourier transforms in the 
impact parameter space of the light-cone wave functions 
and are given by 

f i=n d 2 k- i=n -> 

1p n (x 1 ,p 1 ,X 2 ,P2, ••■) = / Y[ 77T^V GX P(*X!^) 

J i=l * 71 ' i=l 

x 4> n (x 1 ,k 1 ,x 2 ,k 2 ,..)(2ir) 2 6(^2k l ). (13) 




FIG. 2: Geometry of two hard collisions in impact parameter 
picture. 

Thus the 2GPD based description of the four —J- four 
processes is equivalent to the representation for the cross 
section corresponding to the simple geometrical picture, 
but instead of a triple integral we now have an integral 
over one momentum A. The 2GPD defined in Eq. [5] 
is useful for calculation of many different processes. At 
the same time the knowledge of the full double GPD is 



5 



necessary for complete description of events with a double 
jet trigger since the pedestal strongly depends on the 
impact parameter B [4|. 

Let us stress that this picture is a natural generaliza- 
tion of the correspondence between momentum represen- 
tation and geometric picture for a conventional case of 
two — > two collisions. Indeed in this case it is easy to see 
that the cross section in the momentum representation 

a 2 = j f( XuP *)f( X2l p2)^di (14) 

has a simple geometric representation 

f - da h 

<r 2 = d 2 p 1 d 2 Bf(x 1 ,p 1 ,p 2 )f{x 2 ,B - pi,p 2 )-^-<ri, 

J15) 

where f(x,p,p 2 ) = ip + (x, p,p 2 )ip(x, p,p 2 ) and ip(x,p,p 2 ) 
is the Fourier transform of the light-cone wave function 
defined above. 

Let us now summarize our results. We have argued 
that there exists the kinematical domain where the four 
— > four hard parton collisions form the dominant mech- 
anism of four-jet production. In this region we calcu- 
lated the cross section, see Eqs. [2JH] and found that it 
can be expressed through new 2GPDS (see Eq. [5j, ex- 
pressed through the light-cone wave functions of the col- 
liding hadrons. These 2GPDS depend on a transverse 
vector A that measures the transverse distance within 
the parton pairs. (Equivalent expressions for these GPDs 



can be easily given in terms of the operator products.) In 
the impact parameter space we derived the widely used 
intuitive geometric picture. We argued that the observed 
enhancement of a four-jet cross section indicates the pres- 
ence of short-range two-parton correlations in the nucleon 
parton wave function, as determined by the range of in- 
tegral over A. The contribution of perturbative correla- 
tions in the appropriate kinematic domain is suppressed. 
The detailed study of the interplay of the contribution of 
hard/soft correlations will be reported elsewhere [l5j . 

It was argued recently [22[ that the cross section can be 
expressed in terms of two-parton distribution functions. 
Our analysis indicates that a more detailed treatment of 
the QCD evolution effects is necessary. We found that it 
is necessary to introduce the new 2-particle 2 GPDs which 
depend on additional parameter A. The parameter A 
expresses the fact that the difference of the transverse 
components of the parton momenta is not conserved and 
therefore different in \in) and {out\ states in the double 
hard collisions. 
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